Magnetic excitation spectra of NaV205 
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O . 1- Introduction 

O 

It has been recently shown [Q] that sodium vanadium oxide is the second example of 
X ■ the inorganic spin-Peierls compound following CuGeOa 0. Similarly to CuGeOs the 
(73 ■ strong decrease of the susceptibility xiT) in a' — NaV205 is found below T = 34 K as a 
consequence of the opening of a spin-gap. The spin-Peierls transition is one of the most 
interesting phenomena observed in low-dimensional quantum spin systems. It occurs 
in crystals containing linear chains of half-integer spin coupled by antiferromagnetic 
exchange interaction. The spin-Peierls transition is analogue of the Peierls transition 
in metal and results in dimerisation of the lattice induced by magnetic fluctuations. 
The lattice dimerisation and the existence of the spin-gap are also confirmed by x- 
ray, inelastic neutron scattering [Q, NMR |^ and ESR ^ measurements. Using these 
techniques 0, H, |], the spin gap was estimated to be in the range from 85 to 98 K. The 
previously published infrared and Raman ^ spectra show the appearance of 
new modes below the transition temperature, but nothing was said about their origin. 



Based on an early measurement of the NaV205 crystal structure |T0[, it has been 



shown that this oxide, grown in single crystalline form under ambient conditions, has an 



2 



orthorhombic unit cell with parameters a = 1.1318 nm, b = 0.3611 nm, c = 0.4797 nm, 
Z = 2 and the space group P2imn (Dl^). Such a crystalline structure assumes two kinds 
of vanadium chains along the b-axis. One is magnetic V^"*" (S* = 1/2) and the other one 
is a nonmagnetic V^^ (5 = chains). These chains are aligned parallel to each other by 
sharing an edge of adjacent VO5 pyramids, thus forming layers in ab plane. Na atoms 
are situated between layers as intercalates. This magnetic structure is expected to be 
truly one- dimensional since the chains caring the spin are separated from each other by 
nonmagnetic chains along 6-axis. The problem with this structure is missing physical 
argument for such a strong charge modulation. The most recent structural analysis 



shows that all vanadium atoms are equivalent and in mixed valent state. Such a crystal 
structure is also predicted from polarized Raman scattering measurements |l^, |13| . Even 
without charge modulation, due to strong electron correlations, it is possible to have 
formation of antiferromagnetic Heisenberg chains fl^. It is also possible that an induced 
zigzag type of charge ordering along the ladders of V-ions is formed [|l^ . However, in this 
paper we consider only 1-D antiferromagnetic Heisenberg model as an explanation for 
the magnetic behaviour of this compound. This is based on the magnetic susceptibility 
measurements of NaV205 powder |]l[ which shows a broad maximum at 350 K, followed 
by rapid decrease below 34 K. Such a susceptibility behaviour can be well described 
as linear Heisenberg antiferromagnet in the high temperature phase and by Bulaevskii 
theory p[ in dimerized phase, indicating that there is no frustration in this material. 

In this work we present low temperature Raman spectra of NaV205 single crystal. 
The assignation of the observed modes is done on the basis of Monte Carlo calculation 
of excited states of Heisenberg dimerized Hamiltonian, without frustration. We found 
a good agreement between measured and calculated magnetic excitations for 6 = 0.048 
and J = 452 K. The singlet-triplet and singlet-singlet excitations as well as two-magnon 
mode are observed to be at 66, 127 and 132 cm^^ 



2. Experiment 

The single crystals of a' — NaV205 with a size of approximately 1 x 4 x 0.5 mm^ are 



obtained as described in |16]. The spectral lines of an argon- ion laser are used as an 
excitation source. The light beam was focused on the sample surface using cylindrical 
lens with an average power of 100 mW. The scattered light was collected with an 
objective, aperture 1:1.4. The monohromator used was a Jobin Ivon model UIOOO, 
equipped with a Pelletier-effect-cooled RCA 31034 photomultiplier with a conventional 
photon counting system. The temperature was controlled with closed-cycle helium 
cryostat (Leybold). 

Numerical calculation were carried out on a 4-processor Origin 2000 computer, 
courtesy of Institute of Physics Computer Facihties (IPCF). 
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3. Raman spectra 

The polarized Raman spectra measured with 514.5 nm laser line at temperature of 15 K 
are presented in Figure |I[ The new modes found at temperatures below transition 
temperature T = 34 K are denoted in the graph. The lowest excitation is observed 
at 66cm~^, with the largest intensity in ab polarized spectra. This mode is observed 
in bb polarization as well. The continuum of magnetic excitations starts at 132 cm^^ 
which is exactly two times 66 cm~^ energy and is observed only in ab configuration. The 
106 cm~^ mode is found to be largest in bb spectrum. The other modes in the spectra 
are high temperature phonons assigned on the basis of valence shell model calculation 
13| . The broad structure present in aa spectra is assigned as electronic continuum 



which appears due to d — d electronic transitions of V^^ ions [T^ . 

To further analyse the origin of new modes we measured the Raman spectra at 
different temperatures from 15 to 30 K, Figure |^. We found that 66 cm~^ mode exhibits 
the shift towards lower energy of ~ 3 cm~^ by increasing temperature. Besides the 
132 cm~^ mode that shows similar shift, we observed additional structure at 127 cm^^ 
in the T = 15 K spectrum, denoted with an arrow. To resolve these two structures we 
performed higher resolution measurements (resolution< 1.5 cm~^), see inset Figure 
By increasing the temperature this part of the spectrum rapidly looses intensity and the 
difference between these two modes is no longer observable. 

Moreover, the existence of the continuum is confirmed by observing the Fano-type 
asymmetry of the 177 cm~^ phonon with energy overlapping with that of the continuum. 
The asymmetry of this mode is becoming less pronounced by increasing temperature 
and finally becoming complete Lorenzian in the T = 30 K spectrum. 

Other modes in the spectra do not show frequency shift and we assigne them as 
phonons becoming active due to doubling of the unit cell. 



4. Monte Carlo calculation 



Here we present the results of the numerical calculation based on Green's Function 
Monte Carlo technique (GFMC) [0]. We assumed that NaV205 can be modelled as 
1-D Heisenberg antiferomagnet with modulated NN coupling and no frustration. The 
relevance of this model comes from susceptibility measurements |^ were it was shown 
that low temperature behaviour can be well described with the Hamiltonian: 

n = . jY.il + {-ir5%s,+^ (1) 

i 

where Si are spin-1/2 operators. Further, we assumed the periodic boundary conditions 
so that only rings with even numbers of spins were studied. We have calculated energies 
of the ground, the first excited triplet {S = 1) and quintuplet {S = 2) states for 
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rings with = 6, 10, 20, 100, and200 spins. All these energies are calculated exactly 
with no approximation. The resulting triplet and quintuplet gaps. At and Ag, are 
shown in Figure ^ for 6 = 0.05. The extrapolation to — > oo is done using a 
At,q{N) = At,q + {A/N + B/N'^ + C/N^)exp{-N/No) function. Similar simpler 
functions were used in to extrapolate gaps obtained by Lanczos diagonalization for 
systems up to A^ = 28. Here, we have chosen more complicated ansatz to fit points for 
smaller A^ also. The estimated error for the extrapolated values of triplet and quintuplet 
gaps are less then 1%. The size of symbols in figures do not represent this error. Further 
details of extrapolation procedure can be found in pO |. 

The extrapolated results for all 6 considered are given in Figure ^. It is found that 
the singlet-triplet gap can be excellently fitted by the simple formula At/ J = 25'^/^, 
0.01 < 6 < 1. This was also found in the recent study |2l|. In the limit 6^1 
this formula coincides with leading three orders of the perturbation theory, p2|]. The 



6 < 0.01 region was of no concern to us, and we do not know were the crossover to the 
known p3| , ^ critical behaviour At/ J ~ ln6\ occurs. 

The singlet-quintuplet gap can not be fitted in such a simple fashion. We have found 
that it lies above the double singlet-triplet gap value, at the energy Aq = 2AT + e. This 
is consistent with the RPA study p5[ were this mode was identified as an anti-bound 
two-magnon state. The anti-boundedness e decreases towards zero by increasing 6 to 1. 
According to the estimated error of our calculation the existence of a quintuplet state 
at energy differenet then 2 At is established for 5 < 0.1, see inset Figure 0. 



5. Discussion 

To begin discussion, first we compare our GFMC calculation with other numerical 
results. The ground state energy is in excellent agreement with Bethe ansatz solution 
26| in the limit 5—^0 and with the results for isolated dimmers when 6—^1 ||22|| . 



Moreover, our values for the triplet gap coincide precisely with the exact diagonalization 



results ^7j. The most important point, which is one of our main numerical results, is 
exact calculation for the lowest lying quintuplet state, S = 2. As allready mentioned 
above, this mode at energy Ag = 2At + e was identified as an anti-bound two-magnon 



state ||25|. According to that study one should also expect singlet, two-magnon bound 
state, at the energy As = 2At — 2e. Such an excited singlet state can not be calculated 
exactly using GFMC but it can be estimated using the above formula. 

The low temperature susceptibility measurements can be well described with 
Bulaevskii theory, indicating the spin-Peierls transition at 34 K. The fit gave value of 
the exchange integral J and dimerization 5 to be J = 441 K and 6 = 0.048. Using 
slightly different parameters J = 452 K and 6 = 0.048 we calculated singlet-triplet and 
singlet-quintuplet transition energies to be At = 66 cm~^ and Ag = 134.5 cm~^. These 
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values are obtained by a linear interpolation between calculated values for 5 = 0.04 
and 5 = 0.05. The value for singlet-singlet transition = 127 cm~^ is obtained 
using above estimate. This value is in good agreement with exact diagonalization study 
which gives singlet to triplet gap ratio of ~ 1.9 In the Raman spectra the 66 cm~^ 
feature coincides with our numerical result for the triplet gap and it is possible that 
this mode comes from one-magnon scattering process. This value is in agreement with 
other measurements of the spin-gap. Moreover, this mode shifts towards lower energies 
by increasing temperature and by decreasing sodium concentration It was proved 
that sodium deficiency strongly suppresses the SP transition |^9|, which is an additional 
proof for the magnetic origin of this mode. Still, behaviour of this mode in magnetic 
field remains to be studied. Some preliminary results show no change of this mode in 
the magnetic fields up to 7 T [Q. 

The two-magnon mode is found to be at 132 cm~^ which is exactly twice the gap 
energy. It appears as a well-defined peak due to the one- dimensionality of the magnetic 
interaction. Then the structure at 127 cm~^ can be described as singlet-singlet excitation 
with boundness 2e = 5 cm~^, which is in agreement with our calculation. 

It is interesting to compare our results with that of CuGeOa which is the first 
inorganic material exhibiting the spin-Peierls transition. The magnetic mode at 30 cm^^ 
observed in CuGeOs is described as excitation from the ground state to singlet state 
with singlet to triplet gap ratio of 1.76. Such a ratio is a consequence of frustration, 
which is quite large in this material. From our Raman spectra we found that in 
NaV205 this ratio is around 1.95, which is consistent with assumption that frustration 
is negligible in this material. Moreover, the broad magnon band, centred at about twice 
the maximum magnon energy is not observed in our Raman spectra which again proves 



existence of no frustration in NaV205 19 



The appearance of magnetic scattering in ab polarized configuration can not be 
explained using exchange scattering mechanism of Fleury and Loudon |3^, since this 
theory predicts the magnetic scattering in configuration parallel to the chains, ie. in bb 
direction. Our previous study of Raman spectra fl^, using different laser line energies 
shows existence of resonance effects which may influence selection rules for magnetic 
scattering. It is also shown that 66 cm~^ feature, in the spectra measured with 647.1 nm 
laser line consists of two modes separated by ~ 2 cm~^. Still, the mechanism for the 
apearance of the second mode is not yet established. It was proposed that this mode, 
as well as, the 106 cm~^ comes from crystal field excitations of V^"*" ions. 

Finally, we conclude that magnetic excitations in the Raman spectra of NaV205 can 
be well described using GFMC calculation of excited states of dimerized 1-D Heisenberg 
antiferomagnet with J = 452 K, 6 = 0.048 and no frustration. 
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Figure captions 



Figure 1. The polarized Raman spectra of NaV205 at T = 15 K. The modes that 
apear at temperatures below T = 34 K are denoted in the graph. 



Figure 2. The ab polarized Raman spectra as a fimction of temperature. The arrows 
indicate position of magnetic modes, see text. Vertical dashed lines are placed to 
clarify the energy shift. Inset: The 15 K ab Raman spectra measured with increased 
resolution. 



Figure 3. The triplet and quintuplet gaps versus 1/A'', for S = 0.05. The lines are 
fitting functions defined in the text. 



Figure 4. The extrapolated values, N ^ oo, for triplet and quintuplet gaps as 
functions of 6. The full line is a function A/J = 2(5^^^. The dashed line is a function 
A/ J = 4(5^/*, see text. Inset: the gap values with error bars for some 5. 
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